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'We show that the fidehty decay between an initial eigen- 
state state evolved under a unitary chaotic operator and the 
same eigenstate evolved under a perturbed operator satu- 
rates well before the limit, where A'^ is the size of the 
Hilbert space, expected for a generic initial state. We provide 
a theoretical argument and numerical evidence that, for inter- 
mediate perturbation strengths, the saturation level depends 
quadratically on the perturbation strength. 
PACS numbers 

Over the past twenty years different phenomenon 
found in quantum systems that have chaotic classical 
analogs have been suggested as appropriate signatures of 
quantum chaos [1] [2] [3] [4] [5] [6]. Peres [4] conjectured 
that the initial rate and behavior of a system's fidelity de- 
cay due to a small perturbation in the Hamiltonian may 
provide an appropriate signature of quantum chaos. This 
signature provides an analog to the sensitivity to initial 
conditions which characterizes classical chaos but, as a 
consequence of strictly unitary evolution, cannot emerge 
in quantum systems. Recent insights [7] [8] [9] have lead 
to a more detailed understanding of this signature. 

For a unitary map, U , the fidelity compares the evolu- 
tion of an initial state under unperturbed and perturbed 
dynamics. The fidelity is given by 



F{n) = m{u^r{Upur\^^)[^ 



(1) 



where Up = exp{—i5V) is the perturbation operator of 
strength 6, and ipt is the initial state of the system. The 
fidelity decay behavior depends not only on whether the 
map is chaotic but also on the initial state of the system 
and the strength of the perturbation. For chaotic sys- 
tems, the fidelity eventually approaches an asymptotic 
level. Here, we focus on the characteristics of this asymp- 
totic level by studying 
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n — ^oo Tl — ^ 
1 



(2) 



For initial random states the fidelity saturates at 1/N 
[7], as we show below./ However, for eigenstates Foo, 
the saturation level, in much larger and depends sen- 
sitively on the perturbation strength, 5. The study of 
initial eigenstate fidelity decay is of particular interest 
since it is equivalent to the survival probability of a sys- 
tem eigenstate under the influence of a perturbation. 



Below we provide theoretical arguments showing a re- 
gion where Foo depends quadratically on the perturba- 
tion strength. We also test this prediction numerically 
on quantum chaotic maps. 

For chaotic systems and weak perturbation strengths 
the initial fidelity decay is Gaussian, as expected from 
perturbation theory and random matrix theory [10]. For 
stronger perturbations the fidelity decay is exponential. 
The regime of exponential decay, known as the Fermi 
Golden Rule (FGR) regime [8], is reached when a, a typ- 
ical off-diagonal matrix element of perturbation Hamil- 
tonian expressed in the ordered eigenbasis of the system 
Hamiltonian, is greater then the average system level 
spacing, A. It has been shown that for some pertur- 
bations the rate of the the exponential decay increases as 

, until saturating at a rate given by the corresponding 
classical system's Lyapunov exponent [11] [8] [12], or the 
bandwidth of the system Hamiltonian [8] . 

Jacquod and coworkers [8] showed that the fidelity 
decay in the FGR regime is related to the local den- 
sity of states (LDOS) for initial eigenstates. We define 
eigenvectors and eigenangles for the unperturbed oper- 
ator U\vj) — exp{—i4>j)\vj), and the perturbed opera- 
tor UpU\v'^) = expi-icj)'^)\v'^.) . The LDOS is the spectral 
density of the original system under transition rules given 
by the perturbation. Hence, it is as a measure of the over- 
lap between perturbed and unperturbed states separated 
by an angle {(/ij - 0'^) 



(3) 



For an initial eigenstate of U, the fidelity decay is the 
Fourier transform of the LDOS 



Fin) = ^i^i ~ (t>'k)expi-i{(t>j - (t>'k)n). 



(4) 



Previous studies suggest that the LDOS of a complex 
system in the regime of strong perturbation is Lorentzian 
[13-15] 



(5) 



with a width of F = 27rCT^/A where cr is a typical off 
diagonal element of the perturbation operator. Thus, 
using the Fourier transform relation, the initial fidelity 
decay is exponential with a rate of F 



F{n) — exp{—Tn). 



(6) 
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r can be rewritten in terms of perturbation strength as 

follows: (7 = \J S^Vj^j^ where is the second moment of 

the matrix elements Vmn- Vmn be estimated by not- 
ing that for chaotic systems the eigenvectors are random, 
and, therefore, = V/N [9] where = N-^ X^^i ^1 
is the variance of the eigenvalues of V. The average level 
spacing. A, is equal to 2'k/N. The rate of the exponential 
decay, F can now be evaluated as F = S^X^. 

We now turn to the study of the saturation level of the 
fidelity decay, Foe- After a certain amount of time we 
would expect the initial state to become evenly spread 
out over a complete set of states. This implies that F^a 
should be of order 1/N. F^o should also be independent 
of the perturbation strength. A weaker perturbation sim- 
ply leads to a longer period of time until the saturation 
level is reached, but the saturation level should remain 
unchanged. For initial states that arc cigcnstatcs of the 
unperturbed system, however, F^o depends on the per- 
turbation strength. Weaker perturbations, even in the 
FGR regime, lead to saturation levels significantly higher 
then l/N. 

Prosen [7] has noted that for initial eigenstates F^ 1 

in the limit of weak perturbation and F^o (4 — /?) /N 
for strong perturbation where /3 = 1 for maps with cir- 
cular orthogonal ensemble (COE) properties and = 2 
for maps with circular unitary ensemble (CUE) proper- 
ties. Here, we provide a theoretical argument and nu- 
merical evidence for a quadratic behavior for F^o of initial 
eigenstates versus perturbation strength for perturbation 
strengths between these two extremes. 

To evaluate the dependence of the saturation level on 
the perturbation strength let us start by expressing the 
fidelity for an initial eigcnstatc. \vm), as 

F{n) = Kt'ml l]a;me— (7) 

where aim = {u'llvm)- The above equation can be sepa- 
rated into a time independent term plus a time dependent 
term 

= |a/m|^ + \aim\^\akmf cos[(0j - 4)n]. (8) 

I Ik 

The time average of the second term goes to zero while 
the first term determines F^o as an inverse participation 
ratio of the overlap between perturbed and unperturbed 
eigenvectors [7]. In other words, the fidelity saturation 
level is simply the sum of the squared elements of the 
LDOS. 

Using equation (8), we recover the ~ l/N saturation 
level of Prosen in the limit of strong perturbation. An ex- 
tremely strong perturbation could cause the initial state 
(though an eigenstate of the system dynamics) to be- 
come evenly spread over all eigenstates of the system, 
such that Kt^ilvm)! — I/VN for all \v'i)\. F^o would then 



be |a/m|'' — For weaker perturbations the sat- 

uration level will depend on the number of contributing 
eigenvectors, \v'i) and the coefficients aim- This is equiv- 
alent to the width of the LDOS under the particular per- 
turbation. 

Hence, to estimate F^c for intermediate strengths we 
must have an idea of the number of contributing per- 
turbed operator eigenvectors \v'i)\ to the initial eigen- 
state, \vm)- This can be estimated by the width of the 
LDOS, |a;mp. We assume all eigenvectors within the 
width F of the approximate Lorentzian shaped LDOS to 
have equal weight. With this approximation, 

Foo a 1/(F7V) = l/iS^J^N). (9) 

Thus, we expect a quadratic dependence of F^o on the 
perturbation strength in the FGR regime until the satu- 
ration level reaches 0(1/A^). 

A similar analysis for an initial random state, = 
Sm^ml^m)' shows that Fqo — 1/-^ for all perturbation 
strengths. For random states the fidelity can be written 
as 

Fin) = I ^c:;,c,ayaLe-»(^;-<^-)p. (10) 

mlj 

Once again, the right hand side can be divided to a time 
dependent term 

mlj m'l'j' 

(11) 

which vanishes under time average, and time independent 
term 

= jT^EK-I'I^'-I' (12) 

mlj 

where, in the above equations, |cjp ~ 1/7V for a random 
state. For any non-zero perturbation strength, the time 
average of time dependent term will go to zero. The time 
independent term is easily seen to be approximately 1/7V 
in the limits of weak and strong perturbation. For in- 
termediate perturbation strengths we can estimate the 
contribution of the time independent term by analyzing 
the LDOS. Again, we approximate the Lorentzian LDOS 
with a rectangle of width F and height 1/F. Contribu- 
tions to the sum will be non-zero only if the jth and mth 
eigenvectors are a distance of less than F/2 from the Ith 
perturbed eigenvector. Hence, for each of the A'' values 
of I there will be F terms |a/j p and F terms |a/mP each 
of magnitude l/F. The value of the time independent 
term is thus 1/A''. 

The above predictions were first tested on random cir- 
cular unitary ensemble (CUE) maps. Random matrix 
theory predicts the behavior of the fidelity decay in both 
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the Gaussian [10] and FGR [16] perturbation strength 
regimes. The use of a random matrix as the evolution 
operator to study dynamical aspects of quantum chaos 
has been done in [9]. 

We assume that our system is composed of a collection 
of two-level subsystems or qubits. The perturbation used 
is a z-rotation of all of these qubits through an angle S 

rig 

U.^He-^''^'^/' (13) 

where nq — log2N is the number of qubits in the system. 
In the context of quantum information processing, this 
perturbation corresponds to an error in the phase of all 
the quantum bits in a quantum information processor. 
We note that this perturbation also arises in quantum 
control studies as a model of coherent far-field errors [17]. 
For this perturbation, CUE maps exhibit exponential fi- 
delity decay and a Lorentzian shaped LDOS [9] as shown 
in the insets of figure 1. 

Figure 1 shows Foe versus perturbation strength for 
CUE maps, using initial eigenstates of the CUE matrix. 
We see that below the FGR regime there is very little de- 
cay while in the limit of strong perturbation Fao = 2/N 
as expected for CUE maps. Between these we see a power 
law decrease of Foo with increased perturbation strength. 
Since the LDOS is Lorentzian the discrepancy seen in 
figure 1 must be due to the approximation made by re- 
placing the Lorentzian LDOS with a rectangle of width 
r. The actual slope of the data is between 1.8 and 1.9. 
The data is compared to — Ccue/{^'^^'^N), where 
the proportionality constant, Ccue = 3.6 is chosen to 
best fit the data. 




FIG. 1. Saturation level versus perturbation strength for 
initial eigenstates of a random CUE map of dimensions 256 
(circles), 512 (stars), and 1024 (x). For weak perturbations 
below the FGR regime the fidelity barely decays. In the limit 
of strong perturbation Foe saturates at 2/N (solid line). For 
intermediate values of 5, Foo is well approximated by the 
estimate of equation (9) with the proportionality constant 
Ccue = 3.6. Foo is obtained by averaging over 2000 map 
iterations starting at iteration n = 2000, well after the initial 
exponential decay. This is averaged over all A'^ initial eigen- 
states. The lower inset shows the initial exponential fidelity 
decay of the CUE map with A'' = 1024 averaged over all 1024 
system eigenstates. The fidelity decay is plotted versus S^n 
so that the exponential decay rates overlap and the satura- 
tion level is easily seen. The perturbation strengths used are 
0.1, 0.2, 0.3, 0.4 (top to bottom). The upper inset shows a 
semi-log plot of the local density of states for a CUE map 
perturbed by a collective bit ^-rotation, 5 = 0.1,0.2,0.3 and 
0.4 (bottom to top). The solid line is a Lorentzian of width 
r = l^mn/A with Vmn determined numerically from the CUE 
map. 

A similar analysis was carried out for random circular 
orthogonal ensemble (COE) maps. Random COE matri- 
ces can be created from CUE matrices, COE = CUE * 
transpose{CU E) [3]. Like the CUE maps, the COE 
maps have no classical analog and we introduce them 
here as models for the behavior of quantum chaotic maps 
with COE eigenvector statistics and energy level spac- 
ings. Figure 3 shows Foo versus perturbation strength 
for COE maps. Again, an approximate quadratic re- 
lationship emerges but with a different proportionality 
coefficient, Ccoe = 5.4. 

The difference in proportionality constants is in line 
with the work of Prosen [7] who, using a random 
matrix theory argument, predicts a ratio of 3/2 for 
F^^^ / F^^ in the hmit of strong perturbation. We ob- 
serve that this ratio holds for all perturbation strengths 
in the FGR regime. The calculated numerical aver- 
age of F^^ / F^^ for the three Hilbert space dimen- 
sions explored with perturbations in the FGR regime is 
1.48-3/2. 

We next study Foo for a quantum system with a well 
defined classical analog, the quantum kicked top (QKT) 
[18] [3]. The QKT is an exemplary model of quan- 
tum chaos and has been used in previous studies of fi- 
delity decay [4] [7] [8]. The QKT is a unitary map 
Uqkt = exp{—iTTjy/2)exp{—ikJ'l/j) acting on a Hilbert 
space of dimension TV = 2j -I- 1. J is the angular momen- 
tum operator in the irreducible representation and k is 
the kick strength. A kick strength of A: = 12 is used which 
is well in the chaotic region of the QKT. Since the QKT 
shows anti-unitary symmetry, it is part of the COE class. 
The QKT has COE-like nearest neighbor level spacings 
[3] and eigenvector statistics [19]. The same perturba- 
tion, the collective z-rotation, is used. 
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FIG. 2. Foo versus perturbation strength for initial eigen- 
states of a random COE map (x) and the QKT with k = 
12 (circles) of dimensions 256, 512, and 1024 (from top to 
bottom). For weak perturbations below the FGR regime the 
fidelity barely decays. In the limit of strong perturbation Foo 
saturates at (solid line). For intermediate values of S, Foo 
is well approximated by the estimate of equation (9) with the 
proportionality constant C'coe = 5.4. The numerical value of 
-Foo is determined in the same manner as for the CUE maps. 

It should be noted that the data for the QKT and 
COE maps are very similar. This is expected in that, as 
has been conjectured and demonstrated in a number of 
works, quantum chaotic systems have statistical [2] [19] 
and dynamic features [9] similar to those of the canonical 
random matrix theory ensembles. 

The QKT is a system with a classical analog and has 
symmetries not found in random matrices. It is interest- 
ing to see what effect these symmetries, or invariant sub- 
spaces have on Foo- To do this, Foo is calculated for the 
oe subspace (odd under 180° rotations around the y-axis 
[4]) of the QKT which has dimension N = j. The results 
are shown in figure 3 and again we see that Foo approx- 
imately follows a quadratic decrease with increased per- 
turbation strength. However, while the saturation level 
at the limit of strong perturbation does reach the ex- 
pected 3 /N at the same perturbation strength as for the 
full QKT, the intermediate perturbation strengths lead 
to a saturation level that is higher then for the full QKT. 
The coefficient Coe is significantly higher than that of the 
CUE or COE maps. 
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FIG. 3. Foo versus perturbation strength for initial eigen- 
states of the oe subspace of the QKT with k — 12 (stars) and 
the full QKT with k = 12 (circles) for Hilbert space dimen- 
sions 256, 512 and 1024 (top to bottom). For weak pertur- 
bations below the FGR regime the fidelity barely decays. In 
the limit of strong perturbation Foo saturates at 3/A'^ (solid 
line). For intermediate values of 5 of the full kicked top, Foo 
is well approximated by the estimate of equation (9) with the 
proportionality constant Ccoe ~ 5.4. However, Foo for the 
oe subspace does not match Foo of the full kicked top for these 
perturbation strengths. 

In conclusion, we have given a theoretical argument es- 
timating the saturation level of fidelity decay, Foo , for ini- 
tial states that are eigenstates of the system for interme- 
diate perturbation strengths. Numerical simulations for 
systems with and without classical analogs agree with the 
theoretical predictions. However, the presence of invari- 
ant subspaces appears to influence the saturation level of 
the fidelity decay. 
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MIT Institute. 
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